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Abstract 

We consider the possibility that ultra-high energy cosmic rays originate from the annihilation of 
relic superheavy dark matter. We find that a cross section of (gav) ~ KT 26 cm 2 (M x /10 12 GeV) 3 / 2 
is required to account for the observed rate of super-GZK events if the superheavy dark matter 
follows a Navarro-Frenk- White density profile. This would require extremely large-/ contributions 
to the annihilation cross section. We also calculate the possible signature from annihilation in 
sub-galactic clumps of dark matter and find that the signal from sub-clumps dominates and may 
explain the observed flux with a much smaller cross section than if the superheavy dark matter 
is smoothly distributed. Finally, we discuss the expected anisotropy in the arrival directions of 
the cosmic rays, which is a characteristic signature of this scenario. 
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1 Introduction 

The nature of the dark matter and the 
origin of the ultra-high energy (UHE) cos- 
mic rays are two of the most pressing issues 
in contemporary particle astrophysics. In 
this paper we explore a possible connection 



blasi@arcetri.astro.it 

2 rainer@sask.usask.ca 

3 rocky@rigoletto.fnal.gov, 
Rocky.Kolb@cern.ch 



between these issues: dark matter is a very 
massive relic particle species (wimpzilla) , 
and wimpzilla annihilation products are 
the ultra-high energy cosmic rays. 

The dark matter puzzle results from 
the observation that visible matter can 
account only for a small fraction of the 
matter bound in large-scale structures. 
The evidence for dark matter is supported 
by mass estimates from gravitational lens- 
ing [1], by the peculiar velocities of large 
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scale structures [2], by measurements of 
CMB anisotropy [3] and by measurements 
of the recession velocity of high-redshift 
supernovae [4]. Constraints from big-bang 
nucleosynthesis imply that the bulk of the 
dark matter cannot be baryonic, and most 
of the matter density of the universe must 
arise from particles not accounted for by 
the Standard Model of particle physics [5]. 

The existence of UHE cosmic rays of en- 
ergies above the Greisen-Zatsepin-Kuzmin 
cutoff [6], E GZ K - 5 x 10 19 eV, is a ma- 
jor puzzle because the cosmic microwave 
background constitutes an efficient obsta- 
cle for protons or nuclei of ultra-high en- 
ergies to travel farther than a few dozen 
Mpc [6,7]. This suggests that the observed 
extremely energetic cosmic rays with E > 
Eqzk should originate in our cosmic neigh- 
borhood. Furthermore, the approximately 
isotropic distribution of arrival directions 
makes it difficult to imagine that nearby as- 
trophysical sources are the accelerators of 
the observed UHE cosmic rays (but for the 
opposite point of view, see [8]). 

An interesting possibility is that UHE 
cosmic rays arise from the decay of some 
superheavy particle. This possibility has 
been proposed by Berezinsky, KachelrieB 
and Vilenkin and by Kuzmin and Rubakov 
[9,10], see also [11] for a discussion of this 
in the framework of string/M-theory, and 
[12] for a discussion in the framework of 
topological defects. The superheavy par- 
ticles must have masses M x > 10 12 GeV. 
Although this proposal circumvents the as- 
tronomical problems, there are two issues 
to address: Some cosmological mechanism 
must be found for producing particles of 
such large mass in the necessary abun- 
dance, and the lifetime of this very massive 
state must be in excess of 10 20 yr, if we 
want these particles to be both dark mat- 
ter candidates and sources of UHE cosmic 
rays. 



The simple assumption that dark matter 
is a thermal relic limits the maximum mass 
of the dark matter particle. The largest an- 
nihilation cross section in the early universe 
is expected to be roughly M^ 2 . This im- 
plies that very massive wimps have such a 
small annihilation cross section that their 
present abundance would be too large if 
the wimps are thermal relics. Thus, one 
predicts a maximum mass for a thermal 
WIMP, which turns out to be a few hundred 
TeV. While a thermal origin for wimps is 
the most common assumption, it is not the 
simplest possibility. It has been recently 
pointed out that dark particles might have 
never experienced local chemical equilib- 
rium during the evolution of the universe, 
and that their mass may be in the range 
10 12 to 10 19 GeV, much larger than the 
mass of thermal wimps [13-16] . Since these 
wimps would be much more massive than 
thermal wimps, such superheavy dark par- 
ticles have been called wimpzillas [16]. 

Since wimpzillas are extremely mas- 
sive, the challenge lies in creating very few 
of them. Several wimpzilla scenarios have 
been developed involving production dur- 
ing different stages of the evolution of the 
universe. 

wimpzillas may be created during 
bubble collisions if inflation is completed 
through a first-order phase transition 
[17,18]; at the preheating stage after the 
end of inflation with masses easily up to 
the Grand Unified scale of 10 15 GeV [19] 
or even up to the Planck scale [20-22]; or 
during the reheating stage after inflation 
[15] with masses which may be as large as 
2 x 10 3 times the reheat temperature. 

wimpzillas may also be generated in 
the transition between an inflationary 
and a matter-dominated (or radiation- 
dominated) universe due to the "nonadia- 
batic" expansion of the background space- 
time acting on the vacuum quantum fluc- 
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tuations. This mechanism was studied in 
details in Refs. [13,23,24]. The distinguish- 
ing feature of this mechanism is the ca- 
pability of generating particles with mass 
of the order of the inflaton mass (usually 
much larger than the reheat temperature) 
even when the particles only interact ex- 
tremely weakly (or not at all) with other 
particles, and do not couple to the inflaton. 

The long lifetime required in the decay 
scenario for UHE cosmic rays is problem- 
atic if the UHE cosmic rays originate in 
single-particle decays. The lifetime prob- 
lem can be illustrated in the decays of 
string or Kaluza-Klein dilatons. These 
particles can be decoupled from fermions 
[25], and therefore decay in leading order 
through their dimension-five couplings to 
gauge fields 

Ci* = ~ F»\F, V \ (1) 

Here, / is a mass scale characterizing the 
strength of the coupling, and in Kaluza- 
Klein or string theory, it is of order of the 
reduced Planck mass mp\ = (SttGn)" 1 ^ 2 — 
2.4 x 10 18 GeV. In heterotic string theory, 
e.g., f = m Pl /V2 ~ 1.7 x 10 18 GeV [26]. 
If the number of the vector fields is da, we 
find the lifetime estimate from the dilaton- 
vector coupling of Eq. (1) 

327T/ 2 

d G m\ 

1.9 x 10- 22 s Zl0 12 GeV\ 3 
d G \ rrifj, ) 



\1.7 x 10 18 GeV; ' 

If there are direct decay channels through 
first order couplings, superheavy particles 
decay extremely fast, even if the coupling is 
only of gravitational strength and dimen- 
sionally suppressed. Superheavy relic par- 



ticles with a sufficiently long lifetime there- 
fore require sub-gravitational couplings or 
exponential suppression of the decay mech- 
anism due to wormhole effects [9], instan- 
tons [10], or magic from the brane world 
[27]. 

Motivated by the attractiveness of the 
decay scenario, we investigate the possibil- 
ity that UHE cosmic rays may result from 
annihilation of relic superheavy dark mat- 
ter. 

Annihilation of dark matter in the halo 
has been a subject of much interest, with 
particular emphasis on possible neutralino 
signatures in the cosmic ray flux. A rea- 
sonable annihilation scenario for UHE cos- 
mic rays is the production of two jets, each 
of energy Mx, which then fragment into a 
many-particle final state including leading 
particles of energies comparable to a signif- 
icant fraction of Mx- 

Assuming that the relic superheavy dark 
matter follows a Navarro-Frenk-White 
(NFW) density profile [28], in Sec. 2 we 
calculate the expected spectrum and the 
annihilation cross section required to ac- 
count for the observed super-GZK events. 
We find that the necessary cross section 
exceeds the unitarity bound unless the 
annihilation involves very large angular 
momentum contributions or non-standard 
mechanisms, like, e.g., in monopolonium 
decay. 

In Sec. 3 we consider the contribution 
from annihilation in sub-galactic clumps 
of dark matter. We find that annihilation 
in isothermal sub-galactic clumps may ex- 
plain the UHE cosmic rays without violat- 
ing the unitarity bound. 

In Sec. 4 we discuss the expected 
anisotropy in arrival direction of UHE 
cosmic rays if the primary source is su- 
perheavy relic particle annihilation in sub- 
galactic clumps. 

Sec. 5 contains our conclusions. 
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2 Annihilation in the smooth com- 
ponent 

In calculating the UHE cosmic ray flux 
from a smoothQ superheavy dark matter 
distribution in the galactic halo, we as- 
sume a superheavy X-particle halo density 
spherically symmetric about the galactic 
center, nx(d) = n x (d), where d is the 
distance from the galactic center. We will 
assume nx{d) is given by an NFW profile 
[28] 



n x (d) 



d(d + d s 



(2) 



Navarro's estimate for the fiducial radius 
d s for the Milky Way is of order 25kpc 
[29]. Dehnen and Binney have examined 
a flattened NFW profile special case 
of a whole class of halo models and give 
a value d s = 21.8 kpc [30]. We will use 
d s = 3d Q = 24 kpc in our numerical esti- 
mates, where d Q ~ 8 kpc is the distance of 
the solar system from the galactic center. 

The dimensionless parameter Nq may be 
found by requiring that the total mass of 
the Galaxy is 2 x 1O 12 M , which gives 



iVn 



8 x 10 55 2.8 xl0~ 9 kpc 3 



M 



12 



M- 



12 



cm 



where M 12 = M x /10 12 GeV. 

For simplicity, we will assume that wim- 
pzilla annihilation produces two leading 
jets, each of energy Mx, while decay of a 
WIMPZILLA produces two jets, each of en- 
ergy Mx/2. The energy spectrum of ob- 
served UHE cosmic ray events from anni- 
hilation is 



4 We refer to a superheavy dark matter dis- 
tribution as smooth if it can be described by 
a particle density n x (d) which decreases uni- 
formly with distance from the galactic center. 
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47r |d - d | 
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Here, dJ\f(E, E^)/dE is the fragmenta- 
tion spectrum resulting from a jet of energy 
-Ejet- For comparison, the energy spectrum 
of observed UHE cosmic ray events from 
WIMPZILLA decay is 



T = 2 



dN(E,E iet = M x /2) 



d 3 d 



dE 
n x {d) 



x 



An |d - d P 



|2 • 



■01 

where Tx is the decay width of the WIM- 
PZILLA. 

Extrapolations and approximations to 
the fragmentation function dN{E, £j e t) / dE 
have been reviewed in [31]. Birkel and 
Sarkar applied the Herwig Monte Carlo 
program to calculate the spectrum from 
WIMPZILLA decay [32], and Fodor and 
Katz employed a numerical integration of 
the DGLAP evolution equations [33]. The 
resulting spectrum in the interesting range 
between 10 19 GeV and 10 20 GeV is similar 
to the spectrum from the modified leading 
logarithmic approximation (MLLA) [34], 
which was employed by Berezinsky et al. 
in their proposal of UHE cosmic rays from 
decay of superheavy dark matter decay [9] . 

We also use the MLLA limiting spectrum 
in the results of Fig. 1. The salient features 
of the results can be understood by making 
the simple approximation that most of the 
content of the jet is pions, with a spectrum 
in terms of the usual variable x = Ej ' E- v 
(of course, < x < 1), 



let 



dM{ 



x 



-x- 3 / 2 (l 
dx 16 

15 



x 



— aT 3/2 (x < 1). 
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Using this fragmentation approxima- 
tion, the scaling of the flux with Mx can 
be found to be 



2 (<ja)M x 2 (annihilation) 
'TxMZ 1 (decay). 



oc M]l 2 T x M x l (decay). 



(4) 



The factor of M]l 2 is from the fragmen- 



tation function, and the factors of M x 2 
or M^ 1 arise from n\ or nx, respectively. 
Therefore, for a given J 7 , the necessary 
cross section scales as M^ 2 in the anni- 
hilation scenario and the necessary decay 
1 li 

width scales as M x in the decay scenario. 

Calculating the resulting UHE cosmic 
ray flux in the annihilation model, and 
comparing it to the similar calculation in 
the decay scenario, we obtain the results 
shown in Fig. 1. 

The shape of the spectrum is deter- 
mined by the mass of the wimpzilla and 
the overall normalization can be scaled 
by adjusting (a^v) or T x . Clearly, in or- 
der to produce UHE cosmic rays in excess 
of 10 20 eV, M x cannot be much smaller 
than 10 12 GeV. In order to provide enough 
events to explain the observed UHE cos- 
mic rays, {(Jav) has to be in the range 
l(r 25 cm 2 to l(r 27 cm 2 . 

This is well in excess of the unitarity 
bound to the /-wave reaction cross section 
[35-37]: 



4tt . , . 1.5 x KT 47 
<t,u<— (2Z + 1)~ 



M? 



12 



M 2 v 



,100 km s -1 9 
x(2Z + 1) cm 2 . 



The unitarity bound essentially states 
that the annihilation cross section must 
be smaller than M^ 2 . However, as empha- 
sized by Hui [37], there are several ways to 
evade the bound. The annihilation cross 
section may be larger if there are funda- 
mental length scales in the problem larger 
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Fig. 1. UHE cosmic ray spectra from su- 
perheavy particle annihilation (upper panel) 
or decay (lower panel). For both figures the 
solid lines are for Mx = 10 12 GeV and 
the dashed lines are for Mx = 10 14 GeV. 
For annihilation, the solid curve is for 
{(Jav) = 6 x 10~ 27 cm 2 and the dashed curve 
corresponds to (<tav) = 10~ 25 cm 2 . In the de- 
cay case, the solid curve is for T^ 1 = 10 20 yr 
and the dashed curve is for T^ 1 = 8 x 10 20 yr. 

than M^ 1 . One example of a fundamental 
length scale might be the physical size of 
the particle. In this case, the scale for the 
annihilation cross section could be related 
to the size of the particle. Another possi- 
bility for the scale of the annihilation cross 
section might be the range of the interac- 
tion. If the interaction has a range A, then 
one might imagine that the scale for the 
cross section is A 2 . 

A related issue is the typical energy of the 
annihilation products. In this paper we as- 
sume that annihilation produces two jets, 
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each with energy approximately Mx- It is 
easy to imagine that with the finite-size ef- 
fects discussed above, there is the possibil- 
ity that annihilation will produce many soft 
particles, rather than essentially two parti- 
cles each of energy M X - 

An example that suits our needs is the 
annihilation of a monopole-antimonopole 
pair. Assume the monopole mass is M M . 
In the early universe, the annihilation 
cross section is roughly M^ 2 [38]. It is 
this cross section that would determine 
the present monopole abundance. How- 
ever, monopole-antimonopole annihilation 
in the present universe is another matter 
[39]. In the late universe, monopoles and 
antimonopoles can capture with a large 
cross section and form monopolonium, a 
monopole-antimonopole Rydberg state. It 
may take as long as the age of the uni- 
verse for the state to lose energy by radi- 
ation, becoming smaller and more tightly 
bound. Eventually the state will have a 
size of the inner core of the monopole (for 
a GUT monopole, about a/M M ). Then 
the monopole-antimonopole annihilates in 
a final burst of radiation. 

The annihilation spectrum of monopolo- 
nium contains many particles. There are 
many low-energy photons radiated during 
the long decay period, followed by a rapid 
burst of high-energy jets (E ~ aMu) in 
the final death throes of annihilation. 

The monopolonium example also illus- 
trates that the early-universe annihilation 
cross section may be many orders of mag- 
nitude different than the present annihi- 
lation cross section. However, for mag- 
netic monopoles the Parker bound on their 
abundance imposes a further constraint 
[40] . This arises from the requirement that 
magnetic monopoles must not cause too 
much depletion of the galactic magnetic 
field and seems to rule out a noticeable 
UHE cosmic ray flux from their annihila- 



tion in the halo [41]. 

The monopolonium example is an ex- 
istence proof that the annihilation cross 
section may be many orders of magnitude 
larger than M x 2 , while annihilation prod- 
ucts can have energy of order Mx- While it 
is an existence proof, it also illustrates that 
the combined conditions of a a ^> M^ 2 and 
{E) ~ M x probably requires some unusual 
interactions of the dark matter. Of course 
we know essentially nothing about the in- 
teractions of dark matter, in particular its 
self-interactions, which may not be com- 
municated to the visible sector. 

We regard the requisite size of the an- 
nihilation cross section to be the most 
unattractive feature of our proposal. 

3 Annihilation in the clumped com- 
ponent 

So far we have assumed that the galac- 
tic dark matter is smoothly distributed. In 
this section we will consider the contribu- 
tion from inhomogeneities in the galactic 
distribution. 

We have modeled the smooth compo- 
nent of the dark matter by a NFW profile, 
Eq. (2). In addition to this smooth compo- 
nent, the dark matter may have a clumped 
component, as suggested by A-body sim- 
ulations. Using the results of these simu- 
lations, the number of subclumps of mass 
M c i per unit volume at distance d from the 
galactic center can be written in the form: 




where Rf is the core radius of the subclump 
distribution in the galaxy, typically of or- 
der 10 to 20 kpc. The mass of the halo of the 
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Galaxy is M# . The normalization constant 
n° d can be calculated by requiring that the 
mass in the clumped component is a frac- 
tion £ of the halo mass M H . From simula- 
tions, £ ~ 10%. The power index, a, may 
also be found from simulations, with the 
result a ~ 1.9. It is then easy to find 



n,. 



(2-a)e 



An{Rff\n{R H /Rf)M H r] 



2-a ' 



where 77 is the ratio of the mass of the 
largest subclump to the halo mass, ex- 
pected to be in the range 0.01 to 0.1, and 
R H is the radius of the halo. 

Up to this point, all the discussion is in- 
dependent of the density profile of the in- 
dividual subclumps. 

Let us now assume that the density pro- 
file of an individual subclump is an isother- 
mal sphere, with radial profile 



n x (r) = n H (d) 



r 



(5) 



where r is the radius of the individual sub- 
clump, defined as the radius at which the 
subclump density equals the density of the 
background halo at the distance d from the 
galactic center. 

If the subclump mass is M d , we can write 

ro 

M d = AiiM x n H (d) J drr 2 (r/r )- 2 , 


from which we obtain 
1 1/3 



r = 



M rJ 



AnM x n H (d) 



(6) 



The rate of annihilation in an individual 
subclump at distance d from the galactic 
center is 

K{M cU d)=A<Kn 2 H {d)(a A v)4 



X 



3Rn 



+ 



1 QlX T 

^—n 2 H (d)(a A v)-^, 

where R m in is the minimum radius of a sub- 
clump (the inner radius where the profile is 
flattened by efficient annihilations). 

One way to estimate R m in would assume 
equality of the annihilation time scale 



T A = [n(r){a A v)\ 1 



1 



r y 



n H {d)(a A v) \r J 
with the radial infall time 

Tf f = r /v(r ) 

= [4nGM x n H (d)]- 1/2 , 

see e.g. [42], where such an estimate was 
applied to the central core of the galaxy. 

In the present setting this would yield a 
ratio 



Rn 



To 



= (o- A v) 



1/2 



T ff= T A 



n H {d) 
AnGM x 



1/4 



and an annihilation rate in a subclump of 
mass M d at distance d from the galactic 
center: 



n(M ch d) = -M d (a A v) 1 / 2 (4nG) 1 / 4 



x 



n H 



3/4 



M x J 



Now we can calculate the contribution of 
the clumped component in a similar man- 
ner as Eq. (3): 



dAf(E : M: 



dE 



cfd 



1 



27r|d-d | 



Mmax 



dM cl n cl (d,M cl )n(d,M d ). (9) 



M„ 



7 



As before we will use an NFW profile for 
nn(d) with d s = 24kpc. 

Using Rf = 15kpc, rj = 0.1, and£ = 0.1, 
we find 



clumped 



:2 x 10 



smooth 



M x \ 
10 12 GeW 



1/2 



X 



'3 x 10- 26 cm 2> 



1/2 



(<7AV) 



To understand the scaling of the events 
from subclumps, we can make the simple 
approximation that one is at the galactic 
center (|d Q | = 0). In the limit that the con- 
tributions come from d <C Rf and d <C d s , 
we find 

.Foe (a A v) 1/2 M x \ 

Now, we take the density profile of the 
subclumps according to a NFW profile, 



n x lr 



n H {d) 



r [r + r s y 
r[r + r. 



n 2 ' 



(10) 



where r s is the fiducial radius of the sub- 
clump and again r is the radius of the sub- 
clump where nx(r) = nn{d). We assume 
here that the ratio ro/r s ^> 1 is constant for 
all subclumps. The mass of the subclump 
can be written in terms of r and r s as 

M d = AnM x J drr 2 n x (r) 
o 



~4irM x n H (d)rQ In 



which provides the core radius r for the 
subclump once we know r s . 

The rate of annihilations for one sub- 
clump is given by 

ro 

n(M ch d)=4n(a A v) J drr 2 n\{r) 



l -rn-i-n 

+ (4n/3)Rf nm (a A v)n 2 H (R min ), 



which becomes 

n(M cl ,d)^^(a A v)^n 2 H (d). 



To calculate T for NFW subclumps, we 
follow a procedure similar to Eq. 9, but 
with 1Z(M c i, d) given above. Again express- 
ing the flux in terms of the flux from the 
smooth component, we find 



^clumped 3 
T- ' 



smooth 

independent of M x or (a A v). Again, this 
result is represented in Fig. 2. 



10 3 



isotherj 




NFW - all M v 



0.01 



0.1 

(v A a)/3xlO- 26 cm 2 



Fig. 2. The ratio of events from the subclump 
component to the smooth component assum- 
ing either an isothermal or a NFW profile for 
the subclumps. This assumes the ratio (8) 

The ratio of the fluxes from the sub- 
clump component and the smooth compo- 
nent are given in Fig. 2. Comparing Fig. 1 
and Fig. 2, we see that for M x = 10 12 GeV, 
to normalize the flux to the observed events 
would require a cross section 10 3 times 
smaller than the isotropic case if the wim- 
pzilla subclumps have an NFW profile. 

For NFW subclumps the calculation of 
the flux has a natural cutoff at small r. 
This is not the case for isothermal sub- 
clumps, and the resulting flux could be very 
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different with a different choice for R m in- 
For example, one might estimate R m i n by 
requiring detailed hydrodynamic balance 
between the infalling and the annihilating 
matter in the central core of the dark mat- 
ter clumps. In general we can write for the 
change in density of a radially symmetric 
dark matter clump 

d d 

—n x (r,t) = —(n x (r,t)v(r,t)) (11) 
2 

+-n x (r,t)v(r,t) 
r 

-n 2 x (r,t)(a A v), 

where v(r, t) is the velocity of radial infall. 

Detailed balance at r = fi mm amounts 
to the requirement 



v(R r , 



-Rminllx (Rmin) {& AV) 



(12) 



Such a requirement was applied in [43] to 
estimate the core radii of neutralino stars. 

For a particle of energy E = we find 
from energy conservation or from the quasi- 
stationary Euler equation in [43] 

v 2 (r) = -2U(r), 

where the potential U (r) in the clump can 
be calculated as a solution of the Poisson 
equation: 



U(r) = 



2nGM x 



r 

ro 



(13) 



x 



J dr' nx(r')r'(\r — r'\ — r — r'). 



This yields 

U{R m in) ~ -47rGM x n H (d)rlln 



To 



Rr, 



and amounts to a ratio 

Rm.i 



ro 



balance 



n H {d) 
72nGM x 



1/2 



much smaller than what found in Eq. (8). 
With this estimate the rate of annihilation 
in an individual subclump at distance d 
from the galactic center is 



K(d, M d ) w 80M d 



2nGn H (d) 



1/2 



, (14) 



and the flux is 

dJ\f(E, M x ) 



clumped 




x 4.2 m s 



2 „-i 



This is by about a factor 10 15 too large, 
and can be reconciled with the observed 
UHE cosmic ray spectrum only if the WIM- 
PZILLA contribution to Qdark amounts only 
to a small admixture ~ 10~ 15 of annihilat- 
ing wimpzillas in sub-galactic dark mat- 
ter clumps. 

Note that this result is independent of 
the annihilation cross section (<tav), and 
therefore also complies with annihilation 
cross sections well below the s-wave unitar- 
ity limit. With the estimate (12) a small an- 
nihilation cross section implies a very small 
densely packed core region of the dark mat- 
ter clumps. The high density increases the 
luminosity of the clumps and compensates 
for the small factor (cr A v) in the local an- 
nihilation rate. 

4 Predicted signals 

The annihilation rate is very sensitive 
to the local density. This implies that if 
the UHE cosmic rays originate from dark 
matter annihilation and dark matter in our 
galaxy is clumped, the arrival direction of 
UHE cosmic rays should reflect the dark 
matter distribution. 

The first possibility we considered is that 
the dark matter distribution is smooth and 
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follows a NFW profile. In this case, the 
galactic center should be quite prominent. 
In Fig. 3 we illustrate the expected angu- 
lar dependence of the arrival direction. The 
galactic center is prominent for the decay 
scenario where T oc nx [12] , and even more 
prominent for the annihilation case where 
T oc n x . 




Fig. 3. The angular dependence of events 
originating from annihilation or decay of the 
smooth component of dark matter that follows 
a NFW profile. In this figure, cos 8 = 1 corre- 
sponds to the direction of the galactic center. 

Now if we assume that there are isother- 
mal dark matter subclumps within the 
galaxy, then the density in the subclumps 
would be larger than the ambient back- 
ground dark-matter density, and events 
originating from subclumps will dominate 
the observed signal. 

The dominance of events from sub- 
clumps has two effects. The first effect is 
that a smaller annihilation cross section is 
required to account for the observed UHE 
flux. The second effect is that there will 
be a very large probability of detecting a 
nearby sub clump. 

In Fig. 4 we present a histogram showing 
the number of occurrences of single and 
multiple events in arrival directions within 
a square degree. We generated many real- 
izations of the expected subclump distri- 
bution, and calculated the flux from the 
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Fig. 4. The probability of finding a certain 
number of events in a one square degree area 
of the sky 

subclumps that would result in about 100 
detected events. We see that the average 
probabilities are quite reasonable, with 
most events in square degree areas with 
only one event, and a few pairs and triplets 
of events. However if we examine individual 
realizations containing about 100 events, 
we see that there is a large probability of a 
large number of events from single nearby 
subclumps. For instance, in the single re- 
alizations shown in Fig. 4, the isothermal 
subclumps result in a square degree bin 
with 14 events and a square degree bin 
with 33 events. If more than 100 events are 
observed with full sky coverage, the signa- 
ture of subclumps will be unmistakable. 

5 Conclusion 

The origin of UHE cosmic rays from an- 
nihilation (or decay) of wimpzillas has 
the attractive feature of the simplicity of 
generating ultrahigh energies: simple con- 
version of rest mass energy. 

UHE cosmic rays from decay of a smooth 
relic superheavy dark matter population 
in the halo requires a lifetime of about 
10 20 years or a correspondingly reduced 
contribution of the superheavy relics to 

^dark ■ 
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In this paper we have examined the 
annihilation scenario: We found that uni- 
tarity limits very likely exclude an appre- 
ciable contribution to UHE cosmic rays 
from standard particle annihilation in the 
smooth halo component, but monopolo- 
nium decay is a possible counterexample 
to this conclusion. 

We also found that annihilation of wim- 
pzilla particles in sub-galactic clumps 
does not necessarily violate low-/ unitar- 
ity limits. Like the decay of wimpzillas 
of lifetime r ~ 10 10 years, this would also 

imply ^wimpziUa < ^ ^dark- 

The annihilation proposal results in sev- 
eral striking predictions. While not dis- 
cussed in this paper, the annihilation sce- 
nario suggests that the bulk of the UHE 
events are photons. This prediction, com- 
mon to the decay scenario, seems to be 
already in conflict with a recent analysis of 
the composition of the old Haverah Park 
data, which indicate that no more than 
of events above 4 x 10 19 eV are con- 



sistent with being photons [44]. However, 
this limit is weaker for events with energy 
above 10 20 eV. Upcoming data with the 
large statistics of events required to deter- 
mine the composition beyond any doubt 
will soon be available and will definitely 
test the viability of all the top down mod- 
els of UHECRs. 

The true characteristic signature of 
the annihilation scenario is the expected 
anisotropy in arrival direction. If the 
dark matter is smoothly distributed in 
the galaxy, the galactic center should be 
prominent. If the dark matter is clumped 
on sub-galactic scales, then the subclumps 
should be visible. 

Thus, the annihilation scenario can be 
falsified by complete sky coverage. The 
Pierre Auger Observatory and future space 
based experiments, such as EUSO [45] and 
OWL [46] will be able to see the galactic 



center, and, by covering the southern sky, 
should be able to pick out subclumps if 
they are present. 
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